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Abstract 

A  six-accelerometer  inertial  tracking  package  currently 
in  intensive  use  in  measuring  living  human  and  human  surrogate 
response  to  impact  acceleration  is  presented  and  discussed. 

The  discussion  includes  an  enumeration  of  the  various  require¬ 
ments  imposed  upon  the  package  as  well  os  its  design  and 
fabrication . 

The  on-site  calibration  facility  is  described,  including 
a  discussion  of  the  procedures  for  routine  calibration  of  the 
pockoges.  Accounts  of  the  data  acquisition  link  from  the 
packoges  through  the  sled  borne  amplifiers  to  the  hybrid 
computer  are  also  included. 

Particular  attention  is  devoted  to  the  theoretical  aspects 
of  this  system.  A  statement  of  errors  is  developed  and  compared  to 
the  various  precision  parameters  of  this  system  and  to  a  general 
estimation  of  the  dynamic  response  envelope  demonstrating  the 
overall  feasibility  of  this  approach. 


THE  NAVAL  AEROSPACE  MEDICAL  RESEARCH  LABORATORY 
DETACHMENT,  New  Orleans,  is  currently  engaged  in  a  series 
of  experiments  to  determine  human  dynamic  response  to  impact 
acceleration.  Using  a  sled  driven  by  a  pneumotic  occelerator, 
human  volunteers  and  such  human  surrogates  as  large  primates 
and  anthropomorphic  dummies  ore  subjected  to  short  durotion 
accelerations  approximating  crash  situations  and  their  responses 
are  monitored  by  a  wide  ronge  of  measurement  systems. 

These  measurement  systems  fall  into  one  of  three  main 
categories:  Physiological,  Optical,  and  Inertial.  The  physio¬ 
logical  systems  include  electroencephalography,  electrocardio¬ 
graphy,  ond  electrooculography.  The  optical  systems  include 
high  speed  cinematography  (1)  and  television.  The  inertial 
systems  include  sled  mounted  accelerometers,  subject  mounted 
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biaxial  rate  gyroscopes,  and  the  subject  mounted  subminiature 
inertial  fracking  packages  that  are  the  topic  of  this  report. 

Systems  in  the  first  of  these  categories  provide  time  histories 
of  the  volunteer  subject's  elecfrophysiological  status  while  systems 
in  the  next  two  categories  provide  time  histories  of  the  impact 
itself  as  well  as  the  kinematic  response  of  instrumented  segments 
of  the  subject's  anatomy.  This  information  is  required  for  the 
design  and  validation  of  anthropomorphic  dummies  and  of  vehicular 
impact  protection  systems,  and  In  the  measurement  of  a  crash 
victim's  ability  to  perform  escape  or  fire  control  maneuvers 
immediately  after-the  impact. 

THE  SYSTEM 

Each  subminiature  inertial  tracking  package  consists  in 
brief  of  six  accelerometers  mounted  on  an  aluminum  body  that 
is  in  turn  fixed  to  the  subject  anatomy.  Three  of  these  acceler¬ 
ometers  are  grouped  near  a  single  point  and  provide  essentially, 
the  translational  acceleration  of  that  point.  Two  more  acceler¬ 
ometers  are  located  about  a  second  point  and,  used  in  conjunc¬ 
tion  with  the  first  three  accelerometers  provide  two  of  the  three 
components  of  angular  acceleration.  A  single  accelerometer  at 
still  another  point  then  yields  the  third  component  of  the  angular 
acceleration.  Of  course,  many  refinements  of  this  discussion  are 
necessary  in  order  to  workup  and  implement  this  system,  but  the 
most  concise  description  is  that  of  the  six-accelerometer  system 
in  a  three-two-one  configuration  as  shown  in  Figure  1. 

These  tracking  packages  are  designed  to  meet  a  wide 
range  of  requirements.  These  can  be  grouped  in  three  general 
categories:  The  first  of  these  being  the  requirements  imposed 
by  the  nature  of  human  kinematic  response  to  impact  such  as 
dynamic  range,  duration,  bandpass,  and  dimensionality;  the 
second  category  consists  of  limitations  on  the  size,  mass, rotational 
inertias,  and  mountings  of  the  devices  so  that  the  wearing  of  the 
systems  produces  no  significant  perturbations  upon  the  kinematic 
response.  The  final  category  is  composed  of  those  requirements 
consistent  with  the  efficient  and  economical  operation  and 
maintenance  of  the  devices. 

The  nature  of  human  kinematic  response  to  impact  accelera¬ 
tion  has  been  the  subject  of  a  number  of  earlier  studies:  The  joint 
Army-Navy-Wayne  State  University  investigation  presented  by 
Ewing,  et  al  (2)(3)(4)at  the  12th  and  13th  Stapp  Car  Crash  Confer¬ 
ences  and  reported  in  NAMRL  Monograph  21;  work  done  at 
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Fig.  1  -  The  six-accelerometer  system,  3-2-1  configuration 


Holloman  Air  Force  Base  and  reported  by  Bendixen  (5)  and 
Clarke  (6);  and  work  done  ot  W.  Stote  University  reported  by 
Mertz  and  Patrick  (7). 

In  eoch  of  these  studies  the  motions  of  head  (and  in  the 
work  reported  by  Ewing  et  al,  the  first  thoracic  vertebral  body 
as  well)  were  measured  in  terms  of  the  rigid  body  motion  of  some 
coordinate  system  fixed  in  the  anatomy.  Since  the  impact  vectors 
were  confined  to  the  midsaggital  plane,  the  assumption  of  symmetry 
reduces  these  motions  to  three  degrees  of  freedom  per  anatomical 
body;  two  degrees  of  translation  within  the  midsaggital  plane  and 
a  single  rotation  about  an  axis  perpendicular  to  this  plane.  Al  - 
though  this  earlier  work  is  concerned  with  impact  vectors  situated 
in  the  subject  midsaggitol  plane,  on  examination  of  the  results 
suggests  the  following  generalizations  concerning  human  kinematic 
response  to  impact. 

1 .  The  domain  of  response  begins  ot  the  initiation  of  impoct 
and  lasts  from  200  to  350  milliseconds. 

2.  The  response  hos  significont  frequency  components  from 
d .  c .  to  30  to  50  Hz . 

The  approximation  of  rigid  body  motion  is  retained  in  this 
current  work,  but  a  global  range  of  impact  vectors  is  anticipated 
invalidating  the  assumption  of  planar  motion.  Therefore,  six 
kinematic  parameters  must  be  monitored  for  each  rigid  body. 

The  problems  encountered  in  the  design  and  implementation 
of  the  anatomical  mountings  for  these  packages  are  too  numerous 
and  complex  to  be  addressed  here.  But  the  limitations  imposed  on 
pockage  size,  mass,  and  rotational  inertia  that  develop  from 
these  same  anatomical  considerations  are  the  most  important 
factors  bearing  on  the  package  design. 
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At  the  outset  of  the  design  phase,  a  choice  was  to  be 
made  concerning  the  inertial  transducer  types  to  be  incorporated 
in  the  packages.  A  study  of  commercially  available  instrumentation 
had  narrowed  the  design  options  to  packages  consisting  solely  of 
subminiature  piezoresistive  accelerometers,  or  packages  consisting 
of  a  number  of  rate  gyroscopes  as  well  as  these  piezoresistive 
accelerometers. 

The  masses  of  the  accelerometers  are  virtually  negligible, 
but  the  masses  of  rate  gyroscopes  range  from  about  80  gms  for  the 
Northrop  Corp.  Biaxial  Rate  Transducer  ta  100  gms  for  the  single 
axis  Fairchild  and  U.S.  Time  Devices.  Accommodations  for  the 
flange  mountings  on  these  transducers  could  be  expected  to  add 
another  15  to  25  gms.  The  result  is  that  if  nominal  values  for  head 
mass  and  head  rotational  inertia  aboitf  the  Y  <^<is,  a$  in  nodding, 
are  taken  as  4000  gms  and  2.30  x  10°  gm-cm  respectively, 

(as  in  ref.  (8)  ),  that  while  penalties  af  only  2%  ta  3%  additional 
mass  are  incurred  per  transducer  the  corresponding  increase  in 
rotational  inertia  is  about  10%  when  the  anatomical  mounting  is 
at  the  mouth.  These  additional  increases  in  rotational  inertia 
were  deemed  unacceptable,  expecially  at  the  upper  end  of  the 
range  of  impact  levels  included  in  the  experimental  design. 

The  requirements  consistent  with  the  operation  and  mainte¬ 
nance  of  these  devices  call  for  the  minimization  of  associated  sled 
borne  equipment  such  as  amplifiers  and  the  like  and  the  minimization 
of  the  channels  of  information  transmitted  via  hard-wire  whip  cable 
to  the  data  collection  systems.  These  requirements  also  call  for 
modularization  of  the  separate  transducers  used  to  instrument  an 
anatomical  body.  The  concept  of  subminiature  inertial  tracking 
packages  is  in  fact  a  marriage  of  convenience  in  which  a  number 
of  elements  are  joined  to  produce  a  single  albeit  complex  entity 
with  associated  calibration  procedures  and  parameters. 

The  initial  expense  necessitated  by  this  modular  concept  is 
greater;  entire  modules  must  be  maintained  as  backup  instrumentation 
rather  than  a  few  extra  transducers.  This  expense  is  offset  by  the 
ease  with  which  these  modules  can  be  deployed  and  interchanged 
and  by  the  savings  in  bookkeeping  necessary  to  document  the  impact 
experiments.  Moreover,  the  minimal  physical  dimensions  of  the 
piezoresistive  accelerometers  has  so  complicated  the  task  of  locating 
them,  that  anything  less  than  permanent  mounting  becomes  uneco¬ 
nomical  . 

The  subminiature  inertial  tracking  packages  consist  of  six 
piezoresistive  accelerometers  fixed  to  a  ribbed,  aluminum  ‘tee1 
in  a  rigidly  controlled  geometrical  configuration.  This  'tee1 
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then  mates  mechanically  in  a  precise  and  reproducible  fashion 
with  the  various  onatomical  mountings. 

The  rib,  which  is  necessary  to  raise  the  natural  frequency 
of  flexural  vibration  above  the  bandwidth  of  the  kinematic  response, 
precludes  the  economical  machining  of  these  ‘tees'  from  aluminum 
stock  in  quantity;  instead  these  'tees'  ore  sand-cast.  Then  only  the 
surface  opposite  the  ribs  need  be  machined  to  mate  with  the  ana¬ 
tomical  mountings  and  scribed  for  the  accelerometer  placement. 

The  'tee'  is  then  sent  to  Entran  Devices  where  the  accelerometers, 
are  custom  mounted  on  the  scribed  surface. 

These  accelerometers  are  availoble  off  the  shelf  in  o  variety 
of  acceleration  ranges  and  are  damped  to  minimize  spurious  reson¬ 
ances  providing  a  bandwidth  from  d.c.  to  500  Hz  or  more. 

Since  these  transducers  are  designed  for  a  wide  rar\ge  of 
applications  in  which  the  standards  of  accuracy  and  linearity  are 
not  as  stringent  as  those  necessary  for  the  tracking  packages,  a 
preselection  process  identifies  the  devices  best  suited  for  this  work  . 
Batches  of  accelerometers  are  sent  to  NAMRL  for  testing  on  the 
in-house  calibration  facility.  Those  meeting  lab  standards  are 
earmarked  for  inclusion  in  the  tracking  packages  and  the  entire 
batch  is  returned  to  the  monufocturer . 

The  earmarked  transducers  are  then  epoxied  to  the  machined 
'tee'  surface  to  within  two  thousandths  of  an  inch  of  the  scribed 
position.  The  temperature  compensation  modules  are  attached 
at  convenient  sites  along  the  rib.  Finally,  the  leads  from  all  six 
accelerometers  are  dressed  out  three  feet  from  the  assembly  ond 
terminated  on  a  single  connector  that  mates  with  the  sled  borne 
ampl  ifiers. 

Various  stages  of  this  work-up  are  shown  in  Figures  2  and  3; 
the  rough  castings,  the  machined  surface  with  and  without  trans¬ 
ducers,  and  the  final  assembly  complete  with  photographic  targets. 
Figure  4  shows  an  instrumented  human  volunteer. 

SYSTEM  CALIBRATION 

These  inertial  tracking  packoges  are  calibrated  at  regular 
intervals  on  the  same  facility  used  to  preselect  the  accelerometers. 
This  calibration  facility,  diagrammed  in  Figure  5,  consists  of  an 
Inland  Controls  800  seriesVate  table  ir.lerfaced  with  the 
EAI  Pacer  600  hybrid  computer  used  in  the  data  acquisition 
system  (DAS). 

This  rate  table  has  a  twenty-four  inch  diameter  table  top 
and  develops  angular  velocities  of  up  to  ten  thousand  degrees 


Hg.  2  *  Tee-rough  listings  and  the  machined  surface 


l  »g.  }  -  Tee-with  transducers  and  the  linal  assembly 


Fig.  4  -  Instrumented  human  volunteer 
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per  second;  two  small  turntables  mounted  on  the  rate  table  are 
used  to  orient  the  tracking  packages  through  a  wide  range  of 
alignments  relative  to  the  table's  centripetal  acceleration. 

The  calibration  procedures  call  for  the  package  to  be 
mounted  on  one  of  the  turntables  and  the  requisite  orientation 
information  to  be  entered  into  the  DAS.  The  DAS  then  powers 
up  the  rate  table  and  takes  readings  at  regularly  spaced  levels 
of  centripetal  acceleration  throughout  the  nominal  ranges  of  the 
package  instrumentation.  This  operation  is  repeated  for  a  number 
of  different  package  orientations. 

After  six  or  more  such  data  sets  have  been  collected,  the 
voltage  output  of  each  transducer  is  compared  to  the  inertial 
accelerations  at  its  inertial  mass  to  obtain  a  best  fit  for  voltage 
bias  and  the  components  of  the  sensitive  axis.  If  any  of  the  acceler¬ 
ometers  deviate  beyond  certain  limits  from  this  best  fit,  the  package 
is  retired  until  that  unit  is  replaced.  Results  based  on  a  number  of 
these  calibrations  are  shown  in  Figures  6,  7,  and  8. 

The  rate  table-computer  system  is  itself  periodically  calibrated. 
A  device  built  into  the  table  develops  an  electrical  pulse  for  each 
revolution.  These  pulses  are  shaped  and  applied  to  a  Hewlett  Pack¬ 
ard  5330  precision  time  frequency  interval  counter  to  assure  that 
the  angular  velocity  of  the  table  is  exactly  that  commanded  by  the 
DAS.  For  acceptable  accelerometers,  i.e.,  those  meeting  certain 
specified  figures  of  merit,  the  final  output  of  the  calibration  pro¬ 
cedure  consists  of  the  following  parameters:  sensitivity,  offset, 
direction  cosines  identifying  true  sensitive  axis  direction  relative 
to  nominal  accelerometer  coordinates  and  rms  calibration  error. 

This  data  is  keyed  to  a  tracking  package  code  number  and  is  main¬ 
tained  in  a  disk  file  in  the  DAS  to  be  retrieved  each  time  a  package 
is  used  in  a  given  experiment. 

DATA  ACQUISITION 

The  overall  configuration  of  the  DAS  used  in  conjunction 
with  the  inertial  tracking  packages  is  depicted  in  Figure  9.  The 
subject  borne  transducers  are  connected  directly  to  a  set  of  sled- 
mounted  signal  conditioner/line  driver  amplifiers,  the  outputs  of 
which  are  transmitted  via  the  whip  cable  to  the  analog  console 
which  filters  and  further  amplifies  the  incoming  signals.  This 
console  also  functions  as  a  central  distribution  point  from  which 
the  now  buffered  and  amplified  six  accelerometer  outputs  are 
transmitted  not  only  to  the  A/D  converter  For  digitization  at 
a  rate  of  ?000  samples  per  second  per  channel,  but  to  auxiliary 
recording  devices  as  well . 
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SENSITIVITY  VARIATION  FOR  12  50G  ACC  EL  EROMETERS 
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ACCaEROMETER  9 

Fig.  6  -  Variation  in  accelerometer  sensitivity 
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Fig.  7  -  Accelerometer  alignments  -  non  orthogonality  and  variation 
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LINEARITY  CHARACTERISTICS  FOR  12  50G  ACCELEROMETERS 

BASED  ON  4  CALIBRATION  SETS 


© 
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Fig.  8  -  Accelerometer  linearity 


Fig.  9  -  The  data  acquisition  system 
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The  gain  of  a  typical  data  path  is  continuously  adjustable 
from  unity  to  500  under  computer  control.  Based  upon  transducer 
sensitivity  and  expected  dynamic  range  of  a  given  experiment, 
the  gain  is  selected  so  thot  signal-to-noise  ratio  is  maximized. 

The  calibration  of  the  inertial  packages  has  been  previously 
discussed;  the  data  paths  are  also  exhaustively  calibrated  immediately 
prior  to  eoch  experiment.  When  the  sled  omplifiers  are  commanded 
to  the  "calibrate"  mode,  all  amplifier  inputs  are  switched  from  the 
transducers  to  a  common  bus  that  provides  the  calibration  signal. 

The  calibrotor,  on  command  from  the  computer,  generates  a  staircase 
(10-step)  signol  which  is  amplified  and  transmitted  by  each  data 
path;  the  calibrator  output  is  also  returned  to  the  computer  where 
a  comparison  is  made  between  it  and  each  path  output.  Based  on 
this  comporison,  a  total  path  gain  and  offset  is  computed  for  each 
channel,  based  on  a  least  squares  fit  of  the  average  of  150  samples 
of  each  step.  An  interesting  and  useful  feature  of  the  calibration 
subsystem  is  its  capability  for  dynamic  optimization  of  the  calibra¬ 
tion  range.  The  program,  utilizing  the  transducer  sensitivity  and 
signal  path  gain  stored  on  data  files,  automatically  adjusts  the 
amplitude  of  the  10-step  staircase  to  span  the  expected  excursions 
of  the  transducer  outputs  during  the  experiment.  Once  calibration 
is  completed,  the  amplifiers  are  automaticolly  reconnected  to  the 
transducers. 

Inasmuch  as  each  active  path  is  calibrated  prior  to  each  run, 
goin  variations  about  nominol  settings  are  inconsequential,  since 
an  exact  gain  (within  measurement  error)  is  computed.  The  data 
path  parameter  of  critical  importance  relotive  to  overall  acquisition 
accuracy  is  goin  stability,  since  a  gain  shift  between  calibration 
and  acquisition  activities  (nominally  5-10  minutes)  would  yield 
errors.  The  gain  stability  achieved  by  the  NAMRL  system  is  de¬ 
picted  in  Figure  10,  which  shows  that  the  variability  is  inconsequen¬ 
tial  . 

THEORETICAL  DEVELOPMENT 

The  first  consideration  of  the  theoretical  backgrounds  of 
the  inertial  tracking  packages  must  be  the  relationship  between 
the  voltage  output  of  an  accelerometer  and  its  acceleration 
environment . 

The  ‘first  assumption  is  that  the  voltage  output  is  a  function 
of  the  translational  inertial  acceleration  of  some  point  in  the 
accelerometer  as  referenced  to  the  accelerometer  axes: 


Volt  F(A  ,  A  A  ) 
x  y  z 
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Fig.  1 0  -  Gain  stability  of  the  data  paths 


Taking  an  expansion  around  A  -  0 


0) 


Volf  =  Volf  I  +  S  A  +  S  A  +  S  A 
o  xx  y  y  z  z 

+  SS  A  2  +  SS  A  A 
xx  x  xy  x  y 


+  SS  A  A  +  .  .  . 

XZ  X  z 


(2) 


+  SSS  A  3  +  SSS  A  2  A  +  .  .  . 

xxx  x  xxy  x  y 


If  the  accelerometer  is  linear  within  some  region  bounded 

by  |  A  |  =  ±  A 
—  o 

then 

Volts  =  Volts  |  +SA+SA+SA  =  S  •  A  (3) 
o  xx  y  y  zz—  — 
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where: 

Volts  |  and  S  ,  S  ,  ond  S  are  the  parameters  identified 
o  x  y  z 

in  the  calibration  procedures. 

Body  referenced  inertial  acceleration  of  a  point  fixed  in  o 
moving  rigid  body  can  be  expressed  as: 

A=A+goxP+gox(goxP)  (4) 

— p  ^o  —  —  —  — 

Where  A  is  the  inertial  acceleration  of  the  body  referenced  frame 

CO  is  angulor  velocity  of  the  body 

GO  is  the  angular  acceleration  of  the  body 

P  is  the  fixed  location  of  the  point  of  interest 
The  body  referenced  coordinates  are  right-handed  and 
orthogonal .  Cross  products  and  pseudo  vectors  are  defined 
according  to  the  right  hond  rule. 

Applying  equation  (3)  to  equation  (4)  the  voltage  output 
of  an  accelerometer  whose  inertiol  mass  is  locoted  at  P  is: 

Volf -Volf!  =  S  •  A  +  (P  x  S)., ,  -  Q  (5) 

o  —  — o  —  —  _ 

Where: 


Q  =  S  •  (GO  x  (  GO  x  P )  ) 

=  (P-S)  (w)2  -  (S*  w)  (P-tt  ) 


(6) 


Mcking  up  six  vectors  from  S  and  (P  x  S), 

and  from  A  and  GO 
—  o  — 

(Volf  -  Volf  I  )+Q  =  (S,  (PxS)MA  ,  w)  (7) 

o  —  —  —  o  — 

Finally,  given  six  such  accelerometers,  a  six  vector  may 
be  made  of  quantities  corresponding  to  the  left  side  of  equa¬ 
tion  (7)  and  a  six  by  six  matrix  of  (S  ,  (P  x  S )  )  yielding: 


(A  ,  w) 


I 

s 

-pn 


41 


(P  x  S  ) 

-n  — n 


(/ol HVol tl)  +  Q 
*  on  n 


(8) 


th 


where  n  represents  the  n  accelerometer 
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If  the  motrix  mode  up  of  the  occelerometer  porometers  is  not 
ill-conditioned  then  this  procedure  will  yield  the  body  oriented 
inertiol  occelerotions. 

ERROR  ANALYSIS 


Assume  two  sets  of  inertiol  occelerotions ,o  physicol  set 
ond  a  computational  set  reloted  os  follows  to  the  voltoge  outputs 


Volts  =  (S  ,  P  x  S  )-(A  ,  to)  -  Q  +  Volt|  (9) 

n  — n  — n  — n  — o  —  n  on 


and 


Volts  =  (S  1  ,  P1  x  S1  )  .  (A1  ,  (01  ) 
n  — n  — n  —  n  —  o*~ 

-  Q  +  Volt  |  ^  +  € 

n  o  n  n 


(10) 


The  primed  quantities  in  equotion  (10)  ore  the  computotionol 
counterports  of  the  physicol  quantities  in  equotion  (9).  €  is 
o  measurement  error  in  determining  the  occelerometer  oufput. 

Equoting  the  right  sides  of  equations  (9)  ond  (10)  ond 
reorronging 


(A 


W1)  =  I  (s  1  ,  P  1  X  S  1  )  I"1 

n  '  n  n 


I  (S  ,  P  x  S  )  | .  (A  ,  «  ) 

_  n  n  n  o 

-  (Q  -  Q  ’)  +  (Volt  |  -Volt  |  1  ) 

n  n  on  on 


**■•] 


(id 


If  S^  ond  P  ^  correspond  to  the  positions  ond  orientations  shown 

in  Figure  1  ond  differ  from  their  physicol  counterports  by  AS 

—  n 

ond  AP  then  equation  (11)  becomes 
—  n 

(aJ  ,  w1)-  (Ao,  w  )+  (0,0,0,a:2a,'3  u^1  , 


+  M  .  (A  ,  w  ) 

o  — O  — 


(12) 


validated  3-d  inertial  tracking  package 


913 


+  M  *  (  CO  2 

9  1 


2  2 
,  OL  CO_,CO_  A'  &  ) 

1  '  2  3'  3  '  3  r 


•  (Volf  I  -  Volf  I 


o  n 


on 


where  M  and  M  are  six  by  six  mafrices  made  up  of  error 


a  9 


terms  arisinq  ouf  of  AS  and  A  P  . 

—  n  —  n 

The  lasf  four  ferms  on  the  right  side  of  equation  (12)  are 
error  terms.  The  last  three  of  these  are  similar  in  scope  and 
magnitude  to  those  that  would  arise  in  the  nine  accelerometer 
system  proposed  by  King,  Padgaonkar,  and  Krieger  (  9  );  by 
Padgaonkar,  Krieger,  and  King  (10);  and  by  Alem  (11). 

The  first  of  these  four  error  terms  is  unique  to  six  accelero¬ 
meter  systems.  It  arises  out  of  the  LBe  of  computationally  determined 
angular  velocities  in  calculating  ‘Q  1  which  is  the  centripetal,  or 
gyroscopic,  acceleration  componenf^eing  acquired  by  each 
accel  er^rno*'c,r  . 


CO  '  =  CO  CO  CO  -  CO  ^  CO  ^  +  F 

1  12  3  2  3  1 


•  1  • 

C0„  =  CO  -  CO 
2  2 


3 


+  F 


2 


03) 


where  F^,  F^,  and  F^  represent  the  contributions  of  the  last  three 
error  terms  in  equation  (12)  and  are  functions  of  the  package  tra¬ 
jectory  and  the  accelerometer  calibration  errors. 


Assume  very  small  errors  in  the  computed  angular  velocities: 


co  =  co  +  e 


(H) 


Applying  equotion  (14)  to  equation  (13)  and  rearranging  yields 
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p  —  CO  p  —CO  p 

2  3el  1  e3 


®2e3 


=  F, 


05) 


e3-“2el’"l  e2  "  el  e2 


=  F„ 


Letting  products  of  Hie  small  e  terms  go  to  zero  yields  a  set  of 
three  linear  simultaneous  differential  equations  which  may  now 
be  examined  to  obtain  an  understanding  of  the  nature  and  growth 
of  error  in  the  computed  angular  velocities. 


0 

CO 

3 

W2 

el 

F1 

i 

e2 

-CO 

3 

0 

-W1 

• 

e2 

- 

F2 

(16) 

*3 

-CO 

2 

"“l 

0 

e3 

F3 

Letting 

CO 

00 


n  fio 


dt 


07) 


Where  n  is  a  unit  vector,  q  the  magnitude  of  ^  is  a  scalar 

function  of  time  and  <f>  is  a  time  integral  of  ft  and  will  be 
o  o 

called,  for  convenience,  angular  displacement. 

Applying  (17)  to  (16),  dividing  through  by  ft  ,  and  invoking 
the  chain  rule  yields 


08) 


* 

el 

* 

e2 

+ 

# 

e3 

0 

n3 

n2 

el 

F,/.o 

-n3 

0 

"nl 

• 

e2 

= 

O 

*  CN* 
LL. 

"n2 

"nl 

0 

e3 

O 

CO 

IL. 

Where  *  indicates  differentiation  with  respect  to^owhich 
must  be  a  monotomical ly  increasing  function  of  time. 

The  related  set  of  homogeneous  equations  is: 


-  n , 


09) 


0 
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For  near  zero  the  eigenvalues  .are 


Suggesting  that  e  is  a  nearly  sinusoidal  function  of  angular 
displacement.  In  the  region  where  n^  approaches  1  however; 


=  0 

=  1 

=  -1 


(21) 


By  putting  the  eigenvalue  equation  in  polar  coordinates 
such  that  n^  corresponds  to  the  polar  axis  it  can  be  shown  that 

is  at  a  local  minimum  and  X  at  a  local  maximum. 

Even  $o,e  has  components  that  vary  as  an  exponential 
function  of  tfie  angular  displacement,  but  since  this  quantity 
is  never  much  greater  than  a  quarter  of  a  revolution  and  is,  in 
fact,  generally  much  less  for  biodynamic  impact  acceleration, 
study  of  these  homogeneous  equations  shows  a  worst  case  error 
of  less  than  a  five-fold  increase  of  the  initial  error  in  angular 
velocity.  Since  the  initial  value  of  angular  velocity  can  be  so 
closely  controlled  this  increase  is  negligible. 

The  nonhomogeneous  components  of  the  solution  to  equa¬ 
tion  (18)  for  axes  normal  to  the  X  axis  are,  as  might  be  expected, 
mostly  well  behaved  sinusoids  along  with  the  types  of  error  that 
generally  occur  in  straightforward  integration.  The  worst  case 
is  again,  rotation  about  the  X  axis. 

Invoking  techniques  from  Kaplan  (12),  the  nonhomogene¬ 
ous  components  for  this  case  are: 


dt 


(22) 
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and: 


Exp  (<t>)f  (  (  F2  +  F3)  Exp  (-<*>  ))dl 


(23) 


Exp  (-<*>)  J  (  (  F2  “  F3  )  Exp  ( <t>  )  dl 


where  expression  (22)  represents  the  on-axis  solution  and 
(23)  represents  the  two  aff-axis  solutions.  For  angular  displace¬ 
ments  on  the  order  of  a  quarter  of  a  revolution  expression  (23) 
serves  to  treble  or  quadruple  the  errors  t(^at  could  be  expected 
if  it  were  not  necessary  to  calculate  Q  using  integrations  of 
the  computed  angular  accelerations.  But  if  these  sources  of  error 
in  F  are  minimized,  as  has  been  done  by:  Preselection  of  the 
inertial  transducers;  periodic  on-site  calibration  of  the  packages; 
routine  calibration  of  the  datapaths  prior  to  each  experiment; 
and  on  line  digitization  of  the  analog  signals,  then  even  the 
worst  case  errors  in  expression  (23)  become  acceptable. 


EXPERIMENTAL  SET  UP 


The  -G  impact  experiments  currently  in  progress  at  this 
laboratory  have  been  described  in  a  number  of  reports.  (1 ) ,  (13) 

The  essence  of  these  descriptions  is  that  human  volunteers 
equipped  with  anatomically  mounted  inertial  tracking  packages 
located  at  the  mouth  and  the  posterior  spinaus  process  af  the 
first  thoracic  vertebral  body  are  exposed  ta  short  duration 
accelerations  approximating  a  frontal  (eye  balls  aut)  crash 
impact. 

The  position  and  orientation  of  the  tracking  packages  re¬ 
lative  to  coordinate  systems  located  in  the  bony  anatomy  of  their 
mounting  sites  is  measured  using  radiographic  techniques  developed 
at  this  laboratory.  This  information  is  then  applied  to  the  mathe¬ 
matical  procedures  already  described  so  that  the  final  output  of  the 
computations  is  the  time  history  of  the  motion  af  the  anatomical 
coordinate  system  relative  to  a  coordinate  system  fixed  in  the 
laboratory. 

RESULTS 

The  results  of  a  number  of  these  impact  experiments  have 
been  prepared  to  demonstrate  the  workability  of  the  inertial  tracking 
packages.  Since  many  of  the  experiments  include  subject  mounted 
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rate  gyroscopes  as  redundant  instrumentation,  a  check  is  available 
on  the  calculated  angular  velocity  obout  the  subject  Y  oxis. 

This  axis  runs  approximately  right  to  left  in  the  body  accord¬ 
ing  to  the  definitions  in  reference  (4)  so  that  the  motion  for  the  head 
produced  by  such  an  angular  velocity  would  be  a  nodding  motion. 
The  rotation  produced  by  -G  impact  takes  place  primarily  about 
this  axis. 

Figure  11  compares  time  histories  of  head  angulor  velocities 
about  this  Y  axis  as  obtained  from  the  rate  gyroscopes  (*)  and  the 
inertial  tracking  packages  (*)•  Since  the  rate  gyroscopes  behave 
as  a  second  order  low-pass  filter  with  a  natural  frequency  of 
37  Hz,  and  a  domping  rotio  of  .5,  o  similor  plot  was  prepared, 
figure  12,  in  which  the  output  from  the  inertial  tracking  package 
was  oltered  to  opproximote  the  effect  of  this  filtering. 

Off-axis  components  of  the  rotational  motion  are  addressed 
by  a  second  technique.  The  sensitivities  of  certoin  tronsducers 
were  incremented  or  decremented  by  .5%  prior  to  cal  culation*  of 
the  trajectories  in  order  to  generate  o  worst  case  error  envelope 
about  the  solution  for  a  10 g,  -G  impact. 

Figures  13  through  21  show^he  three  body-oriented  angular 
velocities,  the  lab  oriented  three  translational  accelerations,  and 
the  three  body  oriented  angular  accelerations  as  measured  and  cal¬ 
culated  for  the  head.  Figures  22  through  30  show  the  some  informa¬ 
tion  for  T-l,  the  first  thoracic  vertebral  body. 

The  divergence  of  the  envelope  (+)  and  (*)  from  the  solu¬ 
tion  ( • )  for  the  head  angular  velocities  ot  250  milliseconds  is 
about  5%  of  the  peak  magnitude  of  the  head  angular  velocity. 

This  divergence  is  much  greoter  thon  thot  shown  for  the  motion 
of  T-l,  and,  since  the  head  moves  through  an  angle  of  82  versus 
15  for  T-l  ,  this  observation  beors  out  the  conclusions  concerning 
the  effect  of  ongulor  displacement  brought  out  in  the  theoretical 
development. 


CONCLUSIONS 


A  set  of  six  accelerometers  is  sufficient  to  measure  rigid 
body  motion  for  certain  types  of  trajectories.  The  requirements 
for  these  trajectories  are  that  angular  displacement  not  be 
much  greater  than  a  quarter  of  a  revolution  and  that  the  accelera¬ 
tions  themselves  have  no  sharp  peaks  as  might  be  produced  by 
striking  other  rigid  objects. 

The  first  of  these  requirements  is  necessitated  by  the  feedback 
caused  by  gyroscopic  acceleration  terms,  and  the  second  by  the 


TIME  -  SECONDS 


Fig.  1 1  -  Comparison  of  angular  velocity  measured  by  rate  gyroscope  (*) 
and  that  derived  from  an  inertial  tracking  package  (  ) 


Fig.  12  -  Comparison  of  angular  velocity  measured  by  rate  gyroscope  (*) 
and  that  derived  from  an  inertial  tracking  package  (  )  in  which  variations 
in  bandpass  are  eliminated 


VALIDATED  3-D  INERTIAL  TRACKING  PACKAGE 


919 


Fig.  13  -  Head  X  angular  velocity  (  )  and  envelope  (*),  (+) 
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Fig.  14  -  Head  Y  angular  velocity 
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Fig.  15  -  Head  Z  angular  velocity 
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Fig.  16  -  Head  X  translational  acceleration 
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Fig.  17  -  Head  Y  translational  acceleration 
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Fig  1 8  -  Head  Z  translational  acceleration 
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Fig.  19  -  Head  X  angular  acceleration 


Fig.  20  -  Head  Y  angular  acceleration 
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Fig.  21  -  Head  Z  angular  acceleration 


Fig.  22  -  First  thoracic  vertebral  body  (T-l)  X  angular  velocity 
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Fig.  23  -  T-l  Y  angular  velocity 
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Fig.  24  -  T-l  Z  angular  velocity 
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Fig.  25  -  T-I  X  translational  acceleration 


Fig.  26  -  T-l  Y  translational  acceleration 
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Fig.  27  -  T-l  Z  translational  acceleration 
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Fig.  28  -  T-l  X  angular  acceleration 
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Fig.  29  -  T-l  Y  angular  acceleration 
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Fig.  30  -  T-l  Z  angular  acceleration 
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need  for  favorable  scaling  In  minimizing  Hie  error  terms  in  the 
calculation.  This  minimization  of  error  terms  is  itself  made  much 
mare  critical  by  the  gyroscopic  acceleration  terms  imposing  the 
need  far  an-site  calibration  facilities  as  well  as  calibrating  sensi¬ 
tivities  as  vectors  rather  than  idealizing  them  in  terms  of  orthog¬ 
onal  alignments. 

Human  dynamic  response  ta  impact  acceleration  for  cases 
where  no  direct  contact  with  obstructions  is  permitted  results  in 
trajectories  falling  within  the  herein  specified  limits.  If  the 
problems  of  anatomically  mounting  such  sets  of  six  accelerometers 
can  be  overcome,  i.e.  those  of  relative  motion,  added  mass,  re¬ 
strictive  strapping,  etc., then  the  output  of  these  transducers  will 
yield  the  kinematic  parameters  of  the  body  motion. 
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